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Abstract 

R-parity violation in the supersymmetric standard model can be the source of 
neutrino masses and mixing. We analyze the neutrino mass matrix coming from 
either bilinear or trilinear R-parity violation and its collider signatures, assuming 
that the atmospheric and solar neutrino data are explained by three active neutrino 
oscillations. Taking the gauge mediated supersymmetry breaking mechanism, we 
show that the lightest neutralino decays well inside the detector and the model could 
be tested by observing its branching ratios in the future colliders. In the bilinear 
model where only the small solar neutrino mixing angle can be accommodated, the 
relation, 10 3 BR(i/e ± r T ) ~ BR(i/^ =t r =F ) ~ BR(i/r :l= r =F ), serves as a robust test of the 
model. The large mixing angle solution can be realized in the trilinear model which 
predicts BR(ue ± r^) ~ BR(i//x ± t = f) ~ BR(i/r ± r = F). In either case, the relation, 
BR(ej'j) <C BR([ijj) ~ BR(rjj), should hold to be consistent with the atmospheric 
neutrino and CHOOZ experiments. 







1 Introduction 



In the supersymmetric standard model, the gauge invariance and renormalizablity allow 
lepton and baryon number violation and thus it may cause too a fast proton decay Such 
a problem is usually avoided by introducing a discrete symmetry. Among various possi- 
bilities, the Z 2 R-parity and Z 3 B-parity have been advocated as they can be remnants of 
gauge symmetries in string theory [ffj. Imposing R-parity has been more popular because 
of its simplicity and the possibility of having a natural dark matter candidate. The second 
option of allowing lepton number violation is also of a great interest since it can generate 
neutrino masses and mixing [[| in an economical way to explain the current neutrino data. 
There is a huge (but incomplete) list of literature investigating neutrino properties in this 
framework ||. 

R-parity violation may lead to a distinctive collider signature that the usual lightest 
supersymmetric particle (LSP), which is typically a neutralino, produces clean lepton (or 
baryon) number violating signals through its decay Q. In a model of neutrino masses 
and mixing with R-parity violation, one can have more specific predictions for various 
branching ratios of the LSP decay, as the structure of lepton flavor violating couplings is 
dictated by the pattern of neutrino mixing determined from neutrino experiments H || |7|. 
This provides a unique opportunity to test the model in the future collider experiments. 
A necessary condition is of course that the LSP has a short lifetime to produce a bunch 
of decay signals inside the detector. In the models we will consider, the total LSP decay 
rate is proportional to the (heaviest) neutrino mass and thus the measurement of the LSP 
decay length could also be useful to test the model. 

It is the purpose of this paper to examine the correspondence between neutrino oscil- 
lation parameters and collider signatures charactering specific models of neutrino masses 
and mixing from R-parity violation. For this, we will consider the bilinear and trilinear 
models to see whether they can accommodate the atmospheric || and solar neutrino os- 
cillations || and the constraint coming from the CHOOZ experiment jnj, simultaneously. 
One of our basic assumptions is the universality of soft supersymmetry breaking terms 
at a high scale, which is usually imposed to avoid flavor problems in the supersymmetric 
standard model. This implies that the lepton flavor violation occurs only in the super- 
potential with bilinear and/or trilinear R-parity violating terms and the supersymmetry 
breaking mechanism is flavor-blind. Then, the tree-level neutrino mass is generated by the 
renormalization group evolution which breaks universality between the slepton and Higgs 
soft terms at the weak scale. As a specific scheme, we will consider the mechanism of 
gauge mediated supersymmetry breaking which solves the supersymmetric flavor problem 



in a natural way [11]. A comprehensive analysis of neutrino masses and mixing in this 
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context has been performed in Ref. fT^.Q 

Under such an assumption, the bilinear model can only realize the small mixing angle 
of solar neutrino oscillations while the trilinear model can accommodate the large mixing 
angle as well. In both cases, we will investigate whether the LSP decay length is short 
enough and what are the predictions for LSP decay signals which could test the model in 
the future collider experiments. Here, another assumption we make is that the LSP is a 
neutralino. Let us remark that a similar analysis has been made in Ref. [[?]] considering 
supergravity models with generic bilinear R-parity violating terms. 

This paper is organized as follows. In Sec. 2, we calculate the "effective" trilinear 
R-parity violating couplings, rotating away the mixing mass terms between the ordinary 
particles and superparticles which arise as a consequence of bilinear R-parity violation. 
Those couplings are relevant for the LSP decay. In Sec. 3, we examine the neutrino mass 
matrix which is generated through renormalization group evolution and various (finite) 
one-loop diagrams. From this, we will make a qualitative analysis to examine the sizes 
of various R-parity violating couplings which are required to explain the current neutrino 
oscillation data. In Sec. 4, we will provide a numerical analysis to determine R-parity 
conserving and violating input parameters with which the atmospheric and solar neutrino 
masses and mixing are realized, in the context of gauge mediated supersymmetry breaking 
models. Calculating the corresponding LSP decay rate and branching ratios of various 
modes, we will find how the model can be tested in the collider experiments. Finally, we 
will conclude in Sec. 5. 



2 Effective R-parity violating vertices from bilinear 
terms 

Allowing lepton number violation in the supersymmetric standard model, the superpo- 
tential is composed of the R-parity conserving Wo and violating W\ part; 

Wo = ^H 1 H 2 + h'rL i H 1 E^ + hfQ i H 1 D c i +h 1 lQ i H 1 U^ 

W l = e^L^ + ^LiLjEl + ^UQjDl (1) 

Among soft supersymmetry breaking terms, let us write R-parity violating bilinear terms; 

V soft = BfiHiHz + BiCifjtLiHi + m 2 L . H ^LiH\ + h.c. . (2) 
For a recent detailed analysis, see Ref. fTJ. 
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It is clear that the electroweak symmetry breaking gives rise to nonzero vacuum expecta- 
tion values of sneutrino fields, 2>j, as follows 0; 

(i>i) mhm + Bi€ifj,tp 

fl * = / TT \ ~ 2 V / 

where m\. Hx = m\. Hx + q/i 2 , tp = tan/3 = (H 2 )/(H 1 ) and m\ = m|. + M|c 2/3 /2. 
In general, there are three types of independent R-parity violating bilinear parameters 
such as €i, cii and B^/B, which give rise to the mixing between the ordinary particles and 
superparticles. That is, neutrinos and neutralinos, charged leptons and charginos, neutral 
Higgs bosons and sneutrinos, as well as, charged Higgs bosons and charged sleptons have 
mixing mass terms which are determined by the above R-parity violating parameters. 
The mixing between neutrinos and neutralinos particularly serves as the origin of the 
tree-level neutrino masses which will be discussed later. Note that the above quantities 
have to be very small to account for tiny neutrino masses. While the effect of such small 
parameters on the particle and sparticle mass spectra (apart from the neutrino sector) 
are negligible, they induce small but important R-parity violating vertices between the 
particles and sparticles, which make the LSP destabilized and generate one-loop neutrino 
masses. The derivation of the induced R-parity violating couplings has been performed 
in many previous works. The usual approach is to take full diagonalizations of enlarged 
sparticle-particles mass matrices with R-parity violating parts so that the vertices in 
terms of the mass eigenstates are obtained directly. 

In this work, we take an alternative but equivalent approach which is useful when R- 
parity violating parameters are small. It is to rotate away only the small R-parity violating 
(off-diagonal) blocks of the particle-sparticle mass matrices, leaving untouched R-parity 
conserving particle or sparticle masses at the diagonal blocks. In this way, we can draw 
the induced (or "effective") R-parity violating vertices in terms of the electroweak/flavor 
eigenstate basis. A merit of this method is that one can clearly see the vertex structure of 
the induced R-parity violating couplings along with the usual trilinear vertices in Wi of 
Eq. (HD added to the usual R-parity conserving Lagrangian. This is nothing but the usual 
see-saw diagonalization, which we summarize as follows. Let us take sparticle-particle 
mass matrix given by 

M A \ 
At M'J 

with A <C M, M' . Then, the approximate diagonalization (valid up to the second order 
of R-parity violating parameters ~ A/M or ~ A/M') can be done with the help of the 
rotation matrix given by 

i-feet -e \ 
et 1 - |ete ) 
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where can be found by solving the relation, A = M6 — 9M', in the leading order of A. 
The upper and lower diagonal blocks are then shifted as M — > M + (6 A 1 " + A0 t )/2 and 
M' — > M' — (O^A + A^@)/2. Note that the neutrino-neutralino mass matrix has vanishing 
sub-matrix for the neutrinos, M = 0, and the above change in M is just the see-saw 
generation of small neutrino masses. For the other particles/sparticles, such changes can 
be safely neglected. After performing such a rotation, we get the "effective" R-parity 
violating vertices in the electroweak/flavor basis. Then, it is quite straightforward to find 
the corresponding couplings in the mass basis following the usual diagonalization of the 
familiar (R-parity conserving) particle/sparticle mass matrices. 

In this paper, we do not repeat to write the mixing mass terms between sparticles and 
particles. Instead, we will present the rotation matrices 9 in terms of the following three 
bilinear R-parity violating variables; 

€i (or ai) , ^ = Oj-ej, r]i = ai-Bi/B. 

In generic supersymmetry breaking models with non-universality, the above three types 
of parameters are independent. But, in the restrictive models imposing the universality 
condition at the mediation scale of supersymmetry breaking, nonzero values of and 
f]i arise as a consequence of renormalization group evolution and thus only two types of 
parameters are independent. In this paper, we usually take and as independent ones. 

Neutrino-neutralino diagonalization 

Rotating away the neutrino-neutralino mixing mass terms (by 9 N ) can be made by 
the following redefinition of neutrinos and neutralinos: 

where (i/j) and (x°) represent three neutrinos (v e , u T ) and four neutralinos (B, W 3 , H®, H 2 ) 
in the flavor basis, respectively. The rotation elements 0^ are given by 

°ij = & c f c /3 ~~ e ^i3 and (5) 

( n, M z s w M 2 c w M l Mz_ 

[C > } " F N [ c 2 w M 1 + s 2 w M 2 ' c 2 w M 1 + s 2 w M 2 ' Sp pi ' ^ pi ) 

where F N = MiM 2 /(c 2 v M 1 + s^M 2 ) + M\s 2 pj 'pi. Here sw = sin9 w and cw = cos 6^ 
with the weak mixing angle 6w 

Charged lepton/chargino diagonalization 

Defining 9 L and 9 R as the two rotation matrices corresponding to the left-handed 
negatively and positively charged fermions, we have 
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where and e\ denote the left-handed charged leptons and anti-leptons, (x ■ ) = (W , H 1 ) 



and (xt) = (W + ,H 2 ). The rotation elements 0f- R are given by 
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Mi 



(V2, 2s. 



R _ 
My/ . 



-rr&cfcp and 



M -(V2(l-^), 



+ 2 c^) 



w 



w ' F c ■ ll 

and F C = M 2 + M^s 2 p/n. 

Sneutrino/neutral Higgs boson diagonalization 
Denoting the rotation matrix by 9 , we get 



where 



Hi 
\H° 2 J 



\ 



hi + + 0&t 
h° 2 + egu* + 



J 



(7) 



(8) 



eg 
el 



-Oi - 7] iSf3 m A [m^ - m^{m A + M z Sp) - m A M z SpC 2 p}/ 'F t 
+r] i spc (3 m 2 A [m% - ml.{m A + MfcJ) + m\M%(?^\j 'F s 

-Vi44 m A M z[ m li - m A c w}/F s 
+ri i s 3 pC l3 m A M z [m 2 ~ t + m 2 A c w }/F s 



(9) 



with Fs = (to~. — m|)(m|. — m^)(m~. — m^) and m^m^ and m# are the masses 
of pseudo-scalar, light and heavy neutral scalar Higgs bosons, respectively. Note that 
m\ = —Bii/cpSp in our convention. For our calculation, we assume that all the R-parity 
violating parameters are real and so are all #'s. We also note that the presence of the scalar 
fields as well as their complex conjugates in Eq. (§) is due to the electroweak symmetry 
breaking, which is expected to be suppressed by the factor M z /m A . 
Charged slepton/charged Higgs boson diagonalization 



Defining 9° as the rotation matrix, we have 



where 



c 



-a,- 



/ ^ \ 

Hi 
\H 2 J 



2 2 1 
S (3 m A\ m 



( Si - figfTf - 9gH 2 \ 
ef - 9gH{ - 9gH 2 
Hi + egh + Ogef 
V H 2 +9ge l + 9ger J 



(10) 



m ~ m n- 



mfLim 2 Di t/3 



'ID 



5 



eg 


= -m 


SrC 


2/2 2 
nTfl a Win- — Tfl tt 
P A\ Ri H 


- i 




p 2 
777, . urrrU ■ 




(™ 2 H- 


-mljimjj. - 


™L) 




(™ 2 h- -^IJW- 




eg 


= +Vi 




s 2 m 2 A m 2 Di 




+ & 


m1n{m 2 Li - m\- 




(m 2 H . 


-<)« - 


™L) 


(m 2 H - -mjjimjj. - 


-™L) 


eg 


= +Vi 




spc m 2 A m 2 Di 




+ 6 


™ e iV(mh ~ m H 


-) 


(m 2 H . 


-<)W- - 




( m H- - m ln)( m H- ■ 


-™& 



Here, m H - stands for the charged-Higgs boson mass, and m 2 Li , m 2 Ri and m 2 Di correspond 
to the LL, RR and LR components of the i-th charged-slpeton mass-squared matrix, 
respectively, and m~ are its eigenvalues. We remark that the appearance of in 0^ 
and 6*^ is due to the rotations which remove the Goldstone modes from the redefined 
neutral and charged slepton fields. 

With the expressions for the rotation matrices in Eqs. (4)— (11), we can obtain the 
effective R-parity violating vertices from the usual R-parity conserving interaction vertices, 
which are relevant to the LSP decays. We list them below by taking only the linear terms 
in #'s which are enough for our purpose. 

X° — v — Z vertices: 



with 



C x o„ z = tYP L Lt Z ^Zl + h.c. 



[cf,4,0,2cf]0 C/3 . 



Ah 



2c 



w 



(12) 



X° — I — W vertices: 



with L 



x°iw 
x°iw 

ij 



r>X°lW 

ij 



p L i4 m + PrR?° IW 



ejW+ + h.c. 



^2 \. C 1 > C 2 V^cf , C% C2 , C±] C.jC/3 



V2 



[0,-v^cf,0,-^]^ 



c /3 



(13) 



X° — v — Hi 2 vertices: 



with L 
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L ij 



'X.i 



PlL^ 1,2 + PrR^ 



V 3 E \\ + h - C - 



—j= [-tw(Qji - e%), (Oji - of 3 ), (twdji - ef 2 ), 0] 
^ [-t w {e% + 0f 4 ), {e% + 9f 4 ),o, {-t w o» + e%)\ 

[-t w 9j 3 , 6j 3 , 0, 0] 

[-twe s j2 , 9j 2 , 0, 0] 



(14) 
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X° — I — Hi 2 vertices: 



with 



T X°IH+ 
L ij 

T X°IH+ 
L ij 



X°lH+2 



IJ 



+ PrRh 



X°lH+ 2 



tjHi 2 + h.c. 



(15) 



^[^2,^2,^4,0] 



P 



ij 

X°IH+ 
ij 



[V2g% + ^f 2 , -h){§l - 9%), h*ef A ] 



-h<&, -gefi 



In Eqs. (12)-(15), the four components inside brackets correspond to the indices i = 
1, • • • ,4 indicating the neutralino states (B, W 3 , H®, H^)-, respectively, as before. Here, 
let us remark that all of the above vertices depend only on the variables & or rji which 
are generated by renormalization group evolution under the universality condition, even 
though the individual elements 9^, 6f 2 , Ofi and 6^ depend on either o L or a^. This fact 
will be important when we study the LSP decay processes. 

Effective LQd, LQu, LLe and uff* vertices: 

In the below, we list the A-like or A'-like couplings which are, however, neither super- 
symmetric nor S77(2)£-symmetric: 



where 



^■LQd — £ ab 



A dl - 

aij 



A 



d3 
aij 



A aij L ai djP L Q bj + A^L^djPLQbj 
+A-tij (djP L L ai Q bj + L ai P L Q bj dfj + A^L ai P L Q bj d* 



(16) 



+ h.c. 



[On, On] hj , 

[#i3> ®i2\ h d j , 



A 
A 



d2 

aij 

dA 



i3> ^3] , 



^[-M£ + <£,V26&] 



where 



'LQu 



>ab 



Kh L aiQ b jPRUj + AaijL^QhjPnUj 



■c 



(17) 



A ul . 

1 '■aij 



"3 ^bj 1 R^ai ~T ±J a jy^ b jJ 



35 nCi 



A 



M- (u?Q bj P R L m + L c ai Q* bj P RUj ) + A^Q aj P*L ai ] + /i.e. 

a; 

0^,0^] //; • A 



-u3 



•i/l 



-■LLe 



£ab 



h l aijL ai ejP L L bj + A^jL^ejPLLbj 



(18) 



+A«. (ejP L L ai L b3 - L c at P L L b ^) + A^-P^e 



~* 
.7 
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+A%i> j L c aj P R L ai + h.c. 
where k% = [fig, fig] ^ e , A% = [fig, fig] ^ e , A®, = [fi£, fig] ^ 

A lt = ~/=[W^il + ®i2, -V2fia] , Aai = -^=[%fi^ - ^i2> _ V^^i] 

In Eqs. (16)-(18), the two components in the brackets correspond to the two states of the 
SU(2)l doublets with indices a,b= 1,2, and L c = (z/, e c ) is defined as an lepton SU(2)l 
doublet while L c = (z>, e c ) is its scalar counterpart. 
Finally, we have 



C vfI = A^ViP R 



\ u i^ - -djd^ - efe) 



+ h.c. (19) 



where A? = y/2g%. 

As one can see, the above vertices are non-supersymmetric and SU(2)l breaking. But, 
among various terms in Eqs. (16) and (18), one can separate out the supersymmetric 
couplings, £j/tj and e%h e ^ leaving all the vertices depending only on £j or rji similarly to the 
vertices in Eqs. (12)-(15). Then, combining those with the couplings in the superpotential 
([]]), we can define the effective supersymmetric couplings as A^- fc = eih d -5jk + A^- fc and 
Ajjfc = €ihj5jk + \ijk- We will see that these couplings determine the quantity & or r]i 
through the renormalization group evolution of the bilinear (soft) terms. 



3 Radiative neutrino mass matrix from R-parity vi- 
olation 

After performing the rotations described in the previous section, the three neutrinos in 
the "weak-basis" get important mass corrections arising from the see-saw mechanism 
associated with the heavy four neutralinos. As is well-known, this gives the "tree-level" 
neutrino matrix of the form; 

M 2 

^ = -^4, (20) 

which makes massive only one neutrino, 1/3, in the direction of £. The other two get masses 
from finite one-loop corrections and thus z/3 is usually the heaviest component. We fix the 
value of m U3 from the atmospheric neutrino data || and thus the overall size of £ = |£| as 
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where Fn is defined in Eq. (5) and its typical value is given by M 2 . Furthermore, among 
three neutrino mixing angles defined by the mixing matrix 

fl \ / cia s 13 \ ( c 12 



U 









c 2 3 

-S23 



■523 
C23 / 








1 






C13/ 



-S12 





S12 
C12 







1/ 



f22) 



with c 



cos 6 ij and s^- = sin 6^-, etc., two angles are almost determined by the tree- level 



mass matrix 



as follows; 



sin 2 26, 



atm 



sin 2 26 chooz 



sin 2 20. 



sin 2 29 




(23) 



The atmospheric neutrino and CHOOZ experiments || [TtJ require sin 2 26 atm rs 1 and 
sin 2 26 c hooz < 0.2. The angle 6*12 can be determined only after including one- loop correc- 
tions and is responsible for the solar neutrino mixing, 9 so i « # 12 , if the neutrino mass 
matrix is to explain the atmospheric and solar neutrino oscillations as will be discussed 
in the next section. 



An important property of the tree-level neutrino mass matrix (|20|) is that it depends 
on the bilinear R-parity violating quantities & = dj — e« which can be re-expressed as 
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Am 2 + ABifii/j ml.^ 



(24) 



where Am 2 = m 2 Hi — m 2 L ., ABi = B — Bi and m|. = m 2 . + M|c 2 /3/2. In other words, 
the nonzero values of £j arise from the mismatch of soft mass parameters for the Higgs 
field H\ and slepton field Li having the same gauge quantum numbers. If one assumes 
the universality condition, one has Am- = ABi — m\ Hi — at the mediation scale 
of supersymmetry breaking, and their nonzero values are generated by Yukawa coupling 
effects through the renormalization group evolution down to the weak scale. 

Under the assumption that the R-parity violating couplings follow the usual hierarchies 
as the quark and lepton Yukawa couplings, that is, = A' 33 and Aj = A,,33 give the 
dominant contributions, the renormalization group equations (RGE) of the bilinear terms 
are given by 



167r 2 4-Am; 
dt 

9 dii . 
167r 2 ei— AB 
dt 



Qh 2 b X b + 2(1 -5 t3 )h 2 T X T 

ei(6h 2 b A b + 2(1 - 5 i3 )h 2 T A T ) + + 2A i /i T A i 

-AB^h + XiK) 



(25) 
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lQ7t 2 jm 2 LiHl = ~(QX[h b X b + 2\ih T X T ) + m 2 LiHl (3h 2 b + (1 + 5 i3 )h 2 r 



+ (6A^ 6 + 2X, i h T )Am 2 - (6A^ 6 A 6 A^ + 2X i h T A T AA i ) 

where t = InQ with the renormalization scale Q, X b = toq^ + m 2 D c + m 2 Hl + A 2 ., X T = 
m L 3 + m £= + + A 2 . Here, A's are the trilinear soft parameters corresponding to the 
h b , h T , X' { and Aj couplings, and finally AA[ = A\ — A b , AAi = Ai — A T . Under the one-step 
approximation, the above RGE can be solved as 

eiAm 2 - m 2 L H = (3A-/i 6 X 6 + Xih T X T ) In — 

e l AB l = (s'KhA'i + ~Xih T Ai) In — ^ (26) 

where M m is the mediation scale of supersymmetry breaking, and mi is a typical stop mass 
scale where we calculate the sneutrino vacuum expectation values. Here, we have defined 
A- = eih+X'i, K = e i (l-5 l3 )/i r +A i , A^I; = e t h b A b + X'^ and A,ij = e t (l-5 i3 )h T A T +XiA t . 
Note that A3 as well as A3 vanish. In gauge mediated supersymmetry breaking models 
where the mediation scale M m is low, the above approximate solution is quite reliable. 

We are ready to discuss the typical sizes of the supersymmetric bilinear, and trilinear 
parameters, and Aj, Aj, (or Aj, Aj) which will be relevant for the study of the LSP decay. 
Assuming that there is no fine cancellation among various terms in Eq. (^) and the term 
with X b gives the largest contribution in (E6l), we obtain 



ZiCp 3 m b X b M„ 



A' ott z v m 



2 h — ■ (27) 

^ 1 ) <rr> ~ <rr> ~ 



"i 



In gauge mediated supersymmetry breaking models [O], the sfermion soft masses are 
determined by gauge-boson/gaugino loop corrections which implies X b /m 2 . ~ 2q;|/q;|. 



Further assuming the supersymmetry breaking scale A s close to M m , we take M m /mi 
47r/a 3 . This gives 

Ci/i 6 orA<~2O&C0. (28) 

When the tree mass matrix gives the atmospheric neutrino mass scale as discussed, we 
get £jCg ~ 10~ 6 and thus 

e< ~ 20£ iC p/h b ~ 2 x 10- 3 /t p (29) 

for Fn = Mz- This shows that the parameters e, and aj can be very large while main- 
taining £j = aj — 6j very small for low tan 

Let us consider another possibility that Aj gives dominant contribution in Eq. (f26|). 
Following the similar steps as above, we obtain 

€ih T or Aj ~ 560 £jCg (30) 
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where we took X T = 3m~. . This implies that the contribution of A» to is comparable to 
that of A- if Aj ~ 30A-. Later, we will see that the large mixing angle explaining the the 
solar neutrino data can be obtained for Ai^ ~ 5A23. 

So far, we neglected the radiative corrections in the determination of vacuum expec- 
tation values of the sneutrino as well as Higgs fields. To obtain reliable minimization 
conditions for the electroweak symmetry breaking, one has to consider the effective scalar 
potential 

V eff = V + V 1 

where V Q is the tree-level potential and V\ = g^-StrA^ 4 (in ^ — |J includes one-loop 
corrections. With R-parity violation, V\ is a function of not only the Higgs fields but also 



the sneutrinos In deriving Eq. (|24]), we neglected V\ and used the tree-level 

minimization conditions for the neutral Higgs and sneutrino fields. Since the nonzero £j's 
are also generated by the renormalization effect, the inclusion of V% in the determination 
of sneutrino vacuum expectation values is crucial, in particular, in the case of a low-scale 
supersymmetry breaking mediation. In gauge mediation models, such one-loop corrections 
can give rise to an order of magnitude change in neutrino mass-squared values which are 
well- measured in the atmospheric and solar neutrino experiments |[L2|| . After including 
such effects, Eq. (^4j) is modified to 

Am? + ABifitp ^ m\ tHl 
< + < + Eg> 

^ Jl^ (31) 



where E Hl = dV x /H{dHi, eJ. } = dVi/HtdU and sg } = dVi/L*dLi @. In our numeri- 
cal calculation in the following section, we include such improvements. 

In order to get the full neutrino mass matrix, one-loop radiative corrections to neutrino 
mass matrix should be included: 

My = Mj ee + m|; op . 

In the below, we will discuss whether the above mass matrix M v can explain the at- 
mospheric and solar neutrino data, simultaneously. The numerical calculation in this 
direction has been performed first in Ref. [0 without including the effect of the one-loop 
effective potential V\ in the context of minimal supergravity models. Inclusion of full one- 
loop corrections has been made in Refs. [0 and ]n| , in gauge mediation and supergravity 
models, respectively. As mentioned, the one-loop mass matrix M loop lifts twofold degen- 
eracy of the tree-level mass matrix and thus all the three neutrinos get masses which are 
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generically hierarchical. It is instructive to compare the tree and loop mass components, 
in order to get an idea about the mass of the neutrino, z/ 2 , which determine the solar 
neutrino mass scale. The largest contribution to M loop usually comes from the one-loop 
diagrams with A^ and Aj (more generally with the induced ones, A^ and Aj) which takes 
the from 

J 8ir 2 mf — mj m~ 

+ (32) 
8tt z m~ x — m~ 2 m~ 2 

where m . and m f . are the sbottom and stau mass eigenvalues, respectively. 

When A^ > Aj so that A^ give dominant contributions to both £j and M loop , one has 

MSl „ 3 }K<^L „ io-3^ (33) 



assuming the relation (|28[) and 2.5// = m~ b = 500 GeV. In this case, the second neutrino 
mass eigenvalue is determined by the sub-leading contribution of A, to either or M loop . 
Thus, we expect m U2 /m U3 < 10~ 3 £g, or equivalently, 

|j < 10-4- (34) 

Note that the solar neutrino experiments require Am^ = 10~ 5 — 10 -10 eV 2 depending on 
the type of solar neutrino oscillation solutions. As the atmospheric neutrino oscillation 
requires Am§ 2 «3x 10~ 3 eV 2 , it would be much easier to get the so-called vacuum oscil- 
lation or the low Am 2 MSW solution. To realize the large mixing MSW (LMA) solution 
which is now strongly favored by the recent SNO data || , a large tan (3 is needed. Such 
a tendency has also been observed by numerical calculations in the context of minimal 



supergravity models |16], |15| . However, under the assumption that A^ > Aj, it is impossible 
to get a large mixing angle for the solar neutrino oscillation due to the CHOOZ constraint. 
This will become clear when we discuss the bilinear model in the next section. 

Let us now consider the opposite case that Aj ^> A^ so that Aj give dominant contri- 
butions to £j and M loop (which is the case when Aj > 30 A ■ as in Eq. (POD), one finds 

~ 8vr 2 &0c 2 M z m\ ~ & 

for fi = 1.5mf = 200 GeV. Note that m^/m^ can be even larger than one and the 
resultant neutrino mass components satisfy M\ x 12 22 > M" 3 as ^3 = 0- Such a case is 
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not favorable as it cannot give a large mixing angle for the the atmospheric neutrino 
oscillation. From the above discussion, we can infer that the atmospheric and LMA solar 
neutrino oscillation can be realized if the couplings Aj and A^ satisfy a relation in-between 
(33) and (35), which will be the case that Aj is moderately larger than A^.f] We will 
analyze the neutrino masses and mixing in such a scheme and its collider signature in the 
following section. 

4 Atmospheric and solar neutrino oscillations and 
LSP decays in GMSB models 

Let us make a numerical analysis to find how the neutrino mass matrix from R-parity 
violation explain both the atmospheric and solar neutrino oscillations and what are the 
corresponding collider signatures coming from the LSP decay. Our discussions are spe- 
cialized in the models with gauge mediated supersymmetry breaking (GMSB) in which 
the universality condition is automatic and thus supersymmetric flavor problems are nat- 
urally avoided. For our discussion, we will take the minimal number of the messenger 
multiplets (5 + 5) and the messenger supersymmetry breaking scale As not too far from 



the mediation scale M m \T%. We will also concentrate on the cases where the LSP (being 
a neutralino) is lighter than the W boson so that only three-body decays are allowed. 

Due to the (effective) R-parity violating couplings introduced in the previous sections, 
the LSP, denoted by x?, decays through the mediation of on/off-shell W, Z gauge bosons, 
Higgses, sleptons and squarks, producing the following three-fermion final states, 

uuu, ulflf, uqq' , ifqq' . 

Here we do not distinguish the neutrino flavors, and the final quark states will be iden- 
tified with jets. The modes, vlflj and l i jj, are of a particular interest since the flavor 
dependence of R-parity violating couplings, which are relevant to the neutrino mixing 
angles, will be encoded in their branching ratios. 

Before performing the numerical analysis, let us make some qualitative discussions on 
the LSP decay. When the LSP is heavier than the W boson, the decay modes Xi ~~ * 
ifW^ will have sizable branching fractions || |6| and measuring them will give a direct 
information on the ratios £f : Q '■ £i from which we can probe the neutrino mixing angles 
2 The LMA solution may also be obtained in the bilinear model if one relaxes the universality condition 

HQ- 
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#23 and 6*i3 through Eq. fl23|) ||. The decay rate of the mode Xi ~ > is given by 

T(W = -^[|Cf| 2 + |C7^| a ]|6| a 4/ a (^/m|;) (36) 

with (7f = -^[N u c? + N 12 (c% - V^cf) + N 13 (c$ - c%) + N u c%] 

Cf = N 12 cf + -^=N u cl 

Here, iVy are the components of the neutralino diagonalization matrix for the LSP and 
I 2 (x) = (1 - xf{l + 2x). Taking £ iCp = 10" 6 , C{ = 1 and m^o = M z , we get r(W) « 
10~ 14 GeV corresponding to the decay length r rs 2 cm. Thus, the measurement of 
BR(eW) : BR(^W) : BR(rW), or : £ 2 : £ 2 , will be certainly feasible in the future 
colliders. Note that the contribution C^m T /F c ~ m T tp/ji can be neglected unless tan/3 
is very large. 

If the LSP is lighter than the W boson, only three-body decay modes are allowed and 
thus the desired decay modes may be too suppressed to be observed. As a comparison 
with the above two body decay, let us consider the process Xi ~~ * kW* — > kf f whose 
decay rate is 



3G|m 5 -o 

T{hWn = -64^ [|C,l|2 + \ C ^Fr\ 2 m 2 4h(rnyM^) (37) 

where I 3 (x) = [12x - 6a; 2 - 2x 3 + 12(1 - x) ln(l - x)}/x 4 . This gives T(kW*) ta 8 x 10~ 17 
GeV for Cf = 1, facp = 10~ 6 and m^o = 50 GeV. If this is the dominant decay channel, 
the total decay length will be r ~ 2.5m making it hard to observe sufficient LSP decay 
signals. However, it will turn out that the dominant LSP decay diagrams involve the 
effective A^ or Aj couplings which make the LSP decay well inside the detector. This can 
be understood from the previous discussions showing that A-, Aj 3> £jCg. Furthermore, 
the corresponding decay modes Xi ~~ > v 33 or u ^t^J are dominated by the diagrams with 
the exchange of the sneutrino or charged slepton (in particular, the right-handed stau) 
which are relatively light. To get an order of magnitude estimation, let us consider the 
decay rate for Xi ~^ vlflj; 

= ^^fl^lVK^m^m^ml) (38) 

i 

where a' = g' 2 '/An and J is a order-one function of the sparticle masses which is normalized 
to be one in the limit of m~ = m% = 0. Taking X i33 = 2 x 10~ 5 , m^o = 50 GeV = 70 

X{ e fc to lM ' Xl e 3 

GeV and J = 1, we get Y{i>iTf) ~ Y{v 3 lfr T ) ~ 10 -14 GeV which corresponds to r ~ 2 cm. 
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As we will see, this is a typical order of magnitude for the total decay rate of the LSP when 
R-parity violation accounts for the atmospheric and solar neutrino masses and mixing. 

Let us now present our numerical results for the two possible schemes of R-parity 
violation: (i) the bilinear model which has only three input parameters €{\ (ii) the trilinear 
model where we introduce five input parameters A^ and \. 

• Bilinear model: a scheme for the SMA solution. 

The bilinear model with the universality condition is known to accommodate only the 



small mixing angle solution (SMA) of solar neutrino oscillations [0, [L5], which is now 
strongly disfavored by the recent SNO data ||. This model is an attractive option as it 
is the minimal R-parity violating model and provides fairly neat correlations between the 
neutrino oscillation parameters and the collider signatures. In this scheme, the effective 
trilinear couplings are given by A^ = ej/if, and A, = €ih T in Eq. (|26| ) and thus both the 
tree and the loop mass matrix takes the form cx e^-. The other flavor dependence 
comes from the hb, h T Yukawa coupling effects, which is weak unless tan (3 is very large. 
Now that the relation Eq. fl33|) is applied here, the determination of the overall size of £ 
in Eq. (H) and two mixing angles #23 and # 13 in Eq. (p3|) holds almost precisely. Thus, 
the atmospheric and CHOOZ neutrino experiments require <C |£ 2 | ~ (£3 1 which can 
be directly translated to the condition |ei| <C |e 2 | ~ \e%\. This leads to the neutrino mass 
matrix structure; Mj^ < M" 2 < M 22 23 3 3 . As a consequence, only a small mixing angle 
for the solar neutrino oscillation can be accounted for in the bilinear model. Indeed, the 
solar mixing angle # 12 is almost fixed by the relation tan6 l 12 ~ tan #13, and thus we get 
the relation sin 2 29 so i « sin 2 26 c h ooz . 

Given Am 2 2 ~ 3 x 1CT 3 eV 2 for the atmospheric neutrino oscillation, Eq. ( |34|) tells us 
that Am 21 < 3 x 10 _9 i| eV 2 . This estimation is by no means exact but can show some 
qualitative features. For instance, it implies that the right value of Am 21 ~ 5 x 10~ 6 
eV 2 is hardly achieved with small tan j3 in the GMSB models under consideration. In our 
numerical calculation, we looked for the SMA solutions varying the parameters tan/3, A5 
and M m as well as two R-parity violating parameters and e 2 defined at the scale M m 
while keeping e 2 = eg. We could find a reasonable parameter space only for tan (3 ~ 10—25, 
limiting ourselves to tan/5 < 25 because the (right-handed) stau becomes the LSP for 
larger tan f3. 

In Tables 1 and 2, we present two typical sets of parameters accommodating the SMA 
solution and the other neutrino data. In the tables, e° denote input values set at the scale 
M m where supersymmetry breaking is mediated. As can be seen, the effective couplings A- 
and Aj are much larger than ^cp and the dominant decay modes are vjj and vll where the 
diagrams with the exchanges of sneutrinos and charged sleptons give main contributions. 
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Setl 


tan/3 = 10 


As = 40 TeV 


M m = 150 TeV 




7.53 x 10- 6 


2.51 x 10~ 4 


2.51 x 10~ 4 


K 


1.40 x 10~ 6 


4.68 x 10~ 5 


4.68 x 10~ 5 


\ 


7.60 x 10- 7 


2.53 x 10" 5 







2.40 x 10~ 7 


8.01 x 10~ 6 


7.86 x 10- 6 


ii 


1.81 x 10~ 7 


6.06 x 10~ 6 


6.63 x 10~ 6 


BR 


e 


n 


T 


vjj 




3.58 x lO" 1 






1.10 x 10~ 6 


1.22 x 10~ 3 


1.18 x 10-3 


vlfr* 


4.56 x 10" 4 


3.12 x lO" 1 


3.27 x 10" 1 




m Y o=49 GeV 


r = 


7.18 x 10" 14 GeV 



(Am§ 1? Amli) = (2.5 x 10" 3 , 6.1 x 10~ 6 ) eV 2 

(sm 2 2e atm , sm 2 2e sol , sin 2 29 chooz ) = (0.99, 0.0018, 0.0017) 



Table 1: A bilinear model with the input parameters, tan/3, As, M m and e°, allowing 
for the SMA solution. The values of e° are set at the mediation scale M m . The effective 
trilinear/bilinear R-parity violating parameters, A^, Xi/^ defined in the text, are shown 
in the upper part. Here, and £j are the tree-level and one-loop improved values, 
respectively. In the lower part are shown the important branching ratios of the LSP with 
mass m^o and its total decay rate T. The three columns correspond to the lepton flavors, 
% — e, n and r, respectively. For the mode ujj, we do not distinguish the neutrino flavors. 
The resulting neutrino oscillation parameters are presented in the last two lines. 
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Set2 


tan = 25 


As = 45 TeV 


M m = 90 TeV 


4 


8.94 x 10~ 7 


2.98 x 10~ 5 


2.98 x 10~ 5 


K 


4.25 x 1(T 7 


1.42 x 10~ 5 


1.42 x 10~ 5 


\ 


2.30 x 10~ 7 


7.67 x 10~ 6 







4.36 x 1(T 7 


1.45 x 10~ 5 


1.45 x 10~ 5 


ii 


4.04 x 10~ 7 


1.35 x 10~ 5 


1.28 x 10~ 5 


BR 


e 


H 


T 


vjj 




1.11 x lO" 1 




Ifjj 


1.55 x 10~ 6 


1.72 x 10~ 3 


1.72 x lO- 3 


vlfr* 


6.22 x 10~ 4 


4.07 x 10" 1 


4.07 x 10" 1 




m^o=58 GeV 


r = 


4.30 x 10" 14 GeV 



(Amii, Am 2 21 ) = (3.0 x 10~ 3 , 4.1 x lO" 6 ) eV 2 

(sin 2 2^ m , sin 2 2^, sin 2 29 chooz ) = (0.99, 0.0019, 0.0017) 



Table 2: Same as in Table 1 but with tan/5 = 25. 
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One of its consequence is that the total decay rate is larger than 10 -14 GeV making the 
decay length smaller than a few cm. We have checked that the total decay rate is in the 
region of T ~ 5 x 10 -14 GeV for the LSP mass m^o = 25 — 80 GeV. That is, the decay 
length is around r ~ 0.4 cm. This has to be contrasted with the supergravity case |7j 
where one typically gets r > 1 cm for a light LSP. It is also worthwhile to note that the 
SMA solution requires ~ 0.03 ~ ei/e 2 and thus the bilinear model typically predicts 
the following relation: 

10 3 BR(i/e ± r T ) ~ BR(z//i ± r T ) « BR(z/r ± r T ) . (39) 

As pointed out in Ref. [0], the modes kjj are of a great interest. Their decay rates 
are dominated by the W exchange diagrams as the contribution of the largest coupling 
A^ 33 giving Xi ~^ ktb is kinematically forbidden and the coupling X' i22 = €ih s gives the 
sub-leading effect compared to £jCg which enters the x~l~W vertices. Therefore, the ratio 

BR(ejj) : BK(//././) : BR(rjj) (40) 

is almost same as the ratio £f : £f : £f to determine 9 a tm and Q c hooz through Eq. (23) as 
in the case of m^,o > M\y. Here, we remark that the branching fraction BR(ejj) is too 
small to be observed in the future linear colliders. Assuming the integrated luminosity 
1000 ftr 1 per year, the branching ratios below 10 -5 would not be feasible 0. However, 
the measurement of BR(ejj) <C BR(/ijj) « BR(rjj) will provide a robust test for the 
bilinear model. 

• Trilinear model: a scheme for the LMA solution. 

Let us now consider a more general situation that both bilinear and trilinear R-parity 
violating terms are present. In this case, it is convenient to rotate away the supersym- 
metric bilinear terms to the trilinear couplings as we defined the effective ones in the 
previous sections. In this way, we are allowed to introduce only five couplings, A£ and A« 
which are related to the third generation quarks and leptons. This would be the simplest 
trilinear model. 

The trilinear model provides a possibility to realize the LMA solution which is most 
favored at present. As discussed before, in order to get the LMA solution, sizable con- 
tributions to 12 22 are needed to enlarge the solar neutrino mixing while keeping the 
hierarchy of Mij < to realize the large atmospheric neutrino mixing. From the nu- 

merical calculation scanning the five trilinear parameters, we find that the LMA solution 
is realized if A lj2 ~ 5A 23 . The conditions, Ai ~ A 2 and A' 2 ~ A3, are needed to get two 
large mixing angles, while the small CHOOZ angle requires X[ < 0.2A 23 . Under such 
conditions, we could not find any restrictions on the GMSB input parameters tan/3, A5 
and M m . 



18 



Set3 


tan P — 5 


As = 40 TeV 


M m = 80 TeV 


K 


1.07 x 1(T 7 


1.07 x 10~ 5 


0.96 x 10~ 5 


\ 


4.07 x 10~ 5 


4.07 x 10~ 5 







2.66 x 10~ 7 


2.99 x 10~ 6 


2.48 x 10~ 6 




2.57 x 10" 7 


2.67 x 10~ 6 


2.66 x lO" 6 


BR 


e 


fi 


r 






7.92 x 10~ 3 




itn 


5.17 x 10~ 7 


6.55 x lO" 5 


4.49 x 10~ 5 




2.35 x lO" 1 


2.35 x lO" 1 


5.22 x 10" 1 




m*o=46 GeV 


r = 


1.20 x 10~ 13 GeV 



(Am 2 31 , Am 2 21 ) = (3.1 x 10~ 3 , 5.0 x lO" 5 ) eV 2 
(sm 2 26 atm , sm 2 26 sol , sin 2 26 chooz ) = (0.99, 0.80, 0.0001) 



Table 3: A trilinear model realizing the LMA solution. Here the couplings A^ and Aj can 
be considered as input parameters defined at the weak scale. The rests are the same as 
in the previous tables. 
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Set4 


tan (3 = 25 


As = 40 TeV 


M m = 80 TeV 


K 


7.45 x 10~ 8 


4.48 x 10- 6 


7.43 x 10- 6 


A, 


1.61 x 10~ 5 


2.82 x 10~ 5 





Si 


4.32 x 10~ 7 


7.09 x 10~ 6 


1.09 x 10~ 5 




1.62 x 10~ 6 


1.02 x 10~ 5 


1.28 x 10- 5 


BR 


e 


n 


T 






1.54 x 10~ 3 






7.08 x 10" 8 


1.91 x 10~ 5 


4.53 x 10~ 5 




1.15 x 10- 1 


3.53 x 10- 1 


5.31 x lO" 1 




m^o=50 GeV 


r = 


5.62 x 10~ 13 GeV 



(Am§ l5 Am 2 21 ) = (3.0 x 10" 3 , 5.0 x 10~ 5 ) eV 2 
(sin 2 2^ m , sin 2 2^ oZ , sin 2 26 chooz ) = (0.91, 0.84, 0.16) 

Table 4: Same as in Table 3 with tan/? = 25. 

In Tables 3 and 4, we present two examples allowing the LMA solution for tan (3 = 5 
and 25, respectively. We fixed M m = 2A5. The total decay rate is found to be in the 
vicinity of T ~ 5 x 10 -13 GeV for all the LSP mass below 80 GeV. Therefore, the decay 
length is of the order 0.4 mm. This enhancement compared to the bilinear case is due 
to the largeness of the Aj couplings which also make the modes vll dominating. As a 
consequence, we infer the following distinct feature of the LMA solution; 

BR(uef) ~ BR(z///f ) ~ BR(urf) (41) 

with the individual branching ratio is larger than 10%. The relation for the atmospheric 
neutrino mixing angle in Eq. (|23| ) is not as exact as in the bilinear case, but it still holds to 
a good approximation as can be seen from the tables. On the other hand, the expression 
for the CHOOZ angle determined from £j's is not applicable any more and we cannot 
draw any conclusive prediction for the value of it. Still, the relation, BR(ejj) BR(/ijj) 
~ BR(rjj), holds to a certain degree, but these branching fractions become as small as 
10~ 5 making it difficult to be measured in the planned colliders. 
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5 Conclusion 



The supersymmetric standard model without R-parity is an attractive framework for the 
neutrino masses and mixing as certain neutrino oscillation parameters can be probed by 
measuring the decay length and various branching fractions of the neutralino LSP in the 
future colliders experiments. Taking two simple models of R-parity violation, the bilinear 
model with three input parameters and the trilinear model with five parameters, we an- 
alyzed the neutrino mass matrix which explains both the atmospheric and solar neutrino 
data and its consequences on collider searches. One of our basic assumptions is the uni- 
versality of soft terms for which we considered gauge mediation models of supersymmetry 
breaking. A notable consequence of such an assumption is that the LSP (lighter than 
the W boson) decays mainly through the (effective) trilinear couplings A^ and Aj and its 
decay length is found to be in the ballpark of r ~ 0.1 cm. 

The observation of the decay modes vlflj and ifjj will be important as they reflect 
the lepton number violating structure of a certain model. The bilinear model which can 
accommodate only the SMA solution of the solar neutrino oscillation predicts the relation, 
10 3 BR(z/e ± r T ) ~ BR(^ ± r T ) ps BR(z/r ± r =F ). The dominant decay modes are found to 
be Xi ~^ l/ A* ±T=F > ^t ± t =f and ujj, which are all of the order 10%. The trilinear model can 
realize the strongly-favored LMA solution which can be tested by observing the dominant 
decay modes, Xi ~^ ^e ± r T , u/i^t^ and z/r ± r =F , satisfying the relation, BR(i/e i r T ) ~ 
BR(z//i ± r T ) ~ BR(z/r ± r T ). In both cases, the relation, BR(ejj) <C BR(yUjj') ~ BR(rj'j'), 
should hold to be consistent with the atmospheric and CHOOZ neutrino data. 

Acknowledgment: EJC and JDP are supported by the KRF grant No. 2001-003- 
D00037. 
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